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ABSTRACT 


This thesis investigates the feasibility of automating 
the design of microprocessor-based digital filters. The 
ability of a prototype design system to successfully produce 
filter realizations is tested. General filter structures and 
programming algorithms are presented. Shortcomings in the 
Current version of the design system are determined. 
Modifications are made as required to support digital filter 
realizations. The feasibility of filter generation is 
demonstrated using realistic examples taken from the 


literature. 
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ieee eo DUCTION 


Advances in integrated circuit and microprocessor 
technology have allowed increasingly diverse applications of 
these devices. The decrease in their size and power 
requirements and the increased market for them has resulted 
in a decrease in costs while yielding increased computing 
power. If this trend continues, an even broader variety of 
applications can be expected. The effects of increased 
complexity and the decrease in component costs will lower 
total hardware costs. However, as shown by Shooman [Ref. 1], 
software costs can be expected to remain high for three 
reasons: new applications will require new programs to be 
written, the replacement of older, existing computers with 
newer versions will require either completely new software 
or, at the very least, modifications to existing software, 
and because programming is highly labor intensive, it will 
be strongly affected by inflation. In fiscal year 1980, 
fifty-seven billion dollars were spent on computer systems. 
Of this amount, thirty-two billion dollars, or fifty-six 
percent, was allocated to software [Ref. 2]. It can be 
concluded that the most costly expenditures associated with 
computer system development and/or modification are software 
related, and this situation can be expected to remain 


unchanged. It is therefore to our advantage to automate 
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software production. The effort, in terms of manpower, 
required to develop a given software implementation, as well 
fameme Monetary expense necessary to support it, can then be 
Gewored £O the research and development of new systems and 
applications. 

iWiewapplication of digital computers to the problems of 
real time control is an increasingly important aspect of 
Sempucer technology, and is representative of the uses which 
Can be expected with technological advances. The design of 
such systems has proven to be a complex and expensive 
Secertaking. The concepts which will reduce this complexity 
and a methodology for the automated production of real time 
control systems have been developed. [Ref. 3] While no 
Semputer CXists which can duplicate the creative process of 
@esion illustrated in Figure 1, there are elements of it 
which can be tasked to the computer for accomplishment. [In 
order to determine what portions of the design process can 
be automated, the methodology used by the designer must be 
studied. The steps he follows usually involve combining 
existing elements or components in such a way that the 
desired system is realized. The components that are used are 
selected from a library of such items which is located 
either in the designer's memory or is contained in a 
Physical listing. Because the components contained in the 
listing are regularly ordered, a computer can be used to 


maintain it, and, when given the specifications of a desired 
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yemecti, It Call select the necessary entries. This portion of 
the design process has been successfully automated [Ref. 4] 
and is illustrated in Figure 2. The proposed design tool 
merges the ideas of automated software generation and 
automated hardware production into one system. The 
description of this system is the subject of chapter two. 

The original intent of the design tool was to automate 
the production of the software and hardware necessary for 
the realization of real time controllers, but if this tool 
is to be a truly useful design aid, it must be applicable to 
a wider variety of problems. The implementation of digital 
filters is suitable for a dedicated microprocessor system, 
meee typical of potential applications problems. 

Peeetal filtering has been the subject of considerable 
research during the past fifteen years. Various implemen- 
tations have been achieved, including hardwired logic, 
special-purpose computers, and general-purpose computers. 
With the development of the microprocessor, and the recent 
increases in wordsize and speed, a new alternative is 
available. The application of the proposed design system 
to the implementation of microprocessor-based digital 
filters is the subject of this thesis. 

The difference equation form of the filtering function is 
Perel-suited LO microprocessor realization. However, the 
algorithm used can be expressed in a variety of ways. 


Chapter three discusses four possible implementations. The 
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mathematical descriptions discussed are chosen for their 
Simplicity and ease of implementation, but should not be 
construed as the only available alternatives. The 
possibilities are virtually unlimited. The filtering 

function can be performed using one of two general methods, 
independent of the algorithm chosen. The first method 

defines the transfer function as the product of the first and 
second order sections. This is the familiar cascade form of 
the digital filter. The alternative method is the parallel 
form, which expresses the transfer function as the sum of 
first and second order sections. [Ref. 5] The methods for the 
general implementation of each of these applications within 
the bounds of the design system are also discussed in chapter 
Geeee., Lhe optimization of the solution is not of concern. 
The best realization, as determined by metrics such as memory 
Use and chip count, is not the current objective of the design 
system. If the implementation produced satisfies the 
requirements of the problem in terms of timing constraints 

ame runetion, the solution is acceptable. 

Mige=apolication of the design system to problems other 
than controller design is important for several reasons. By 
demonstrating the ability of the system to generate success- 
ful realizations for a wide variety of problems, its overall 
utility as a design tool will be proven. Varying the type 
of implementation problem presented to the system extends 


bhe applicability of the design technique and refines our 
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understanding of the system. By testing various applications, 
inadequacies in the current implementation of the system can 
be detected and corrected. The attempt to generate digital 
filters using the design system revealed a number of such 
problems and the search for these inconsistencies became a 
major portion of the thesis. The problems found are presented 
met corresponding solutions in chapter four. 

Chapter five describes the application of the design 
system to realistic problems. The concepts discussed in the 
previous chapters are utilized to implement both the cascade 
and parallel forms of a fourth order digital filter. 

The system described in this thesis represents a useable 
design tool. Unfortunately, the coding that the user is 
required to produce in order to achieve the desired hardware 
and software realization is both awkward to use and tedious 
to generate. The addition of new realization volumes, as 
well as the maintenance of existing ones, is a difficult 
process as well. The potential of the design system is 
fmeretore limited by aur ability to simplify the user input 
specification and provide a means by which the software and 
hardware library may be modified to accommodate any potential 
design problem. Recommendations for these and other 


improvements are contained in the sixth and final chapter. 
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iio ube OVERY LEW 


The current version of the design system has evolved 
from the model proposed by Matelan [Ref. 6] and implemented 
by Ross [Ref. 7]. Each of these efforts used the design of 
real time controllers as the problem model for development 
@f the system. The concepts and terminology introduced by 
Matelan can be found in the current implementation. The 


results of this work are summarized here. 


fee ROBLEM MODEL 

The first aspect of the system to be considered is the 
Beeplem model. In order to produce a successful realization 
Tor any problem, it must first be expressed in a form 
understandable by the design system. 

meeecontingency/Task Pairs 

Pach problem Can be expressed as a collection of 

Ccontingency/task pairs, where a contingency is defined as a 
Beeiecal or relational function of some input variable or 
variables. The associated task is executed when the 
contingency is satisfied. Therefore, the first step in 
developing a realization is to express the problem in terms 


of contingency/task pairs. 
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2. High Level Problem Description 

Implicit in the model of the problem are rigid 
constraints on the testing of each contingency and the time 
allowed for execution of the corresponding task. The 
designer must specify these real time requirements in the 
statement of the problem. Matelan proposed a new high level 
design language called a Control System Design Language, or 
CaP, tor this purpose. The language consists of four 
sections: identification section, environment section, 
contingency list, and procedures section. The identification 
Beetion 1S a record of the user identification of the 
problem. The environment section specifies the interface 
between the device and the process which is to be operated 
Meom.e. Lt defines all input and output variables and their 
Characteristics, as well as those variables internal to the 
Mathematical operation of the device. The contingency list 
is a declaration of those conditions that the device must 
respond to, the associated task that each must execute, and 
the real time constraints imposed upon each pair. The timing 
constraints are determined by the maximum time allowed to 
recognize a contingency and the maximum time available to 
execute the corresponding task. The routines which comprise 
each contingency/task are contained in the procedures 
section. By definition, contingencies are written as 
functions, while tasks are specified as procedures. Each 


@emiains the high level descriptions necessary for the 
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performance of its role in the system being produced. Figure 
3 is an example CSDL listing of a simple filtering problem. 
The identification section is readily found and easily 
understood. The environment section specifies two variables, 
one for input and one for output. Each contains eight 
TTL-compatible bits. The contingency section indicates that 
the function READY is executed every millisecond, and when 
true, the task FILTER is performed. The procedures section 
contains the descriptions of the steps necessary to perform 
Pm@emrest for the contingency READY and the execution of the 
process FILTER. In the case of the contingency READY, an 
external variable, DATA, is read. If its value is one, 
indicating the presence of data for processing, the 
contingency is satisfied and the value of READY is set to 
@meeeand the task FILTER is executed. If, on the other hand, 
DATA is equal to zero, the contingency will not be satisfied 
and READY will remain equal to zero to indicate a false 
Semdation. The task FILTER is not executed under these 
Circumstances. 

iMiceruneuron FILTER contains the description of the 


difference equation 
page= xm) + O0.5y(n-1) - 0.5y(n=2). 


The value of the input is assumed to be in digital form and 
memreadd first. The value of the corresponding filtered 


output is then computed and the result is issued as an output. 


ILS, 





BOENTIFICATION: 
DESIGNER: M, R.~ Heilstedt 
DATE: 19 April 1983 
PROJECT: CSDL Example - Filter 


ENVIRONMENT: 
rHewTs X,o,7Tl 
OUTPUT: Y,8,TTL 


CONTINGENCY LIST: 
wine READY:1MS Do FILTER 


PROCEDURES: 
Contingency READY 
Sense (DATA); 
If DATA=!1 Then READY:= fi; 
Exit READY 


Imaisk FILTER 
Sense (X); 
Y=GEeoeoe 1 Ni J=(0.5*YNe)-; 
Issue (Y); 
YNeC=YN13 
YN1I=Y; 

Exit FILTER 


aaiene 3 ee EeccdD Lo te SDE. bac cane 


20 





EE ———— 


The current version of the design system does not 
support che use of CSDL for problem specifications, requiring 
Mifemene User generate the list of primitives manually. A 
user specification format based on the Ada programming 
language is currently under development. [Ref. 8] In the 
interim, CSDL serves as a useful tool in understanding the 
meemeture Of the problem specification required by the 
design system and the transformation it undergoes in the 


course of generating a design. 


fee RMEDIATE FORM OF THE PROBLEM 

The input specification is translated into an inter- 
femmate representation consisting of a list of primitives. 
This transformation is analogous to the operation performed 
by a compiler on high level languages such as Fortran and 
Pascal. The list of primitives is a complete description of 
the problem, containing all of the information necessary to 
generate a hardware and software implementation of the design 
problem and to verify that all timing requirements have been 
Satisfied with the selected realization volume. 

A second file, called the IJADEFL file, is generated at 
the same time as the primitive list, and contains the list 
of contingency/task pairs and their corresponding timing 
Semstraimts. This data is extracted from the ID table, the 
environment table, the timing data from the contingency 


list, and the design criteria table. The design criteria 





table is generated from the data listed in the design 
criteria section of the CSDL listing. This section was 

added by Ross to provide a method for the designer to specify 
the criteria upon which an acceptable realization could be 
chosen. The criteria are: first realization generated, 
least costly realization, and the realization that uses the 
masc power. The current version of the system only supports 
the first realization generated criterion. The creation of 
mem rADEFIL file and the primitive list is the first action 
taken by the design system in its attempt to generate a 
hardware realization. The operation of the désign system in 


Metmaston to these files is depicted in Figure 4. 


Sees UNCTIONAL ELEMENTS 

iiesprobléem formulation is performed by the user for 
Botm the high level description and intermediate represen- 
Metron. Each of these tasks will eventually be automated 
as well, providing a user friendly interface to the design 
Seeeem. AS stated previously, development of the input 
module is in progress, but until a satisfactory input 
specification to intermediate form translator can be 
developed, the input to the design system must be manually 
compiled. 

1. Optimizer Module 

Piter the input files have been produced, they are 


Passed to the first component in the design system, the 
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optimizer module. It serves two purposes, acting as the main 
Meoeream Controlling the operation of the system, and func- 
meomeme as the input module for the primitive list. Although 
Heopeimplemented in the current version of the system, this 
module will have the capability of comparing the multiple 
realizations generated by the system and selecting that one 
which best satisfies the chosen metric as specified by the 
Seem che design criteria section of the input listing. 
feeetunctional Mapper 

The Functional Mapper is the first module called by 
mienwoOpctimizer, and is used to ensure that each primitive 
contained in the intermediate problem specification can be 
M@eatized with the current realization library. This is 
accomplished using two separate operations. The realization 
volume index is first searched for the primitive name as 
eeen by the current line of the intermediate list. If the 
primitive is found, the specifications associated with the 
Meamitive are compared to those contained in the realization 
tibrary. The mapping is considered successful only if ail of 
the primitives can be realized and their criteria satisfied. 
If a given primitive cannot be mapped to a realization, or 
meats selection criteria fail to fit those of the realization, 
a new realization library will be searched. Failure to 
satisfactorily realize a primitive in any library results in 


an unsuccessful mapping. 
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oe taming Analyzer 


The second module executed is the Timing Analyzer. 
It generates the monitor necessary for controlling the 
operation of the device under development, ordering the 
execution of the realization contingency tests and task 
executions such that the timing constraints will be satisfied. 
In order to make such an assurance, the timing analyzer 
assumes that all contingencies are true. This requires that 
all of the tasks must be executed by the monitor and there- 
fore defines the worst case situation. The resulting premise 
is that if the worst case can satisfy the timing requirements, 
all cases satisfy them as well. The timing analyzer determines 
the length of time required to execute all of the code con- 
tained in each of the contingency/task pairs and compares it 
to the timing constraints listed in the Applications Timing 
Teme as generated from the IADEFL file. If all of the 
memanes Constraints are met, the realization is successful. 
imeemot, the realization fails. The output of the Timing 
Pmelyzer 1s used to generate monitor primitives for success- 
ful realizations. These primitives determine the sequence of 
execution for the contingency/task pairs, and are added to 
the list of primitives derived from the high level description 
of the problem. 

An important result of the research conducted by Ross 
was that the design system is capable of automatically 


producing multi-processor realizations. This is done in the 
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Timing Analyzer. If a single processor realization is 
impossible, the Timing Analyzer partitions the Applications 
Timing Table into two separate lists. Each of these contains 
the timing information corresponding to complete sets of 
contingency/task pairs. In other words, the partitioning of 
the problem results in groupings of these pairs. The 
regularly ordered nature of the contingency/task model makes 
this possible. In order to generate separate sets of hard- 
ware, the two timing tables are passed to the Formatter 
module individually. Theoretically, it 1s possible to pro- 
duce a realization in which each contingency/task pair is 
located in its own processor. It is important to note that 
such a system may require shared resources. The current 
version of the design system will only test for a dual 
processor realization when a single one fails. The 
generation of such a realization cannot be forced. 

we Frormatter Module 

The Formatter receives the complete list of 

primitives necessary to produce the output listing for the 
design. The listing consists of the software necessary to 
perform the desired function and the hardware required to 
support it. The hardware and software listings are written 
tO separate files which are in turn copied to an output 
device such as a terminal or line printer. The design 


program is then terminated. 
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D. REALIZATION LIBRARY 

The realization library has thus far been mentioned only 
in describing the operation of the other components of the 
design system, but its importance in the determination of a 
successful design should not be overlooked. Each realization 
library contains volumes of hardware and software primitives 
Dased on individual processor families. The only volume 
currently implemented uses the Intel 8080 microprocessor and 
its various support devices to generate realizations. 

The general format of each line in a volume is the same. 
Fach line is assigned a unique identifying number, and 
contains a maximum of eighty characters. The first set of 
lines within the volume serve as an index to the primitives 
it contains. The lines of the index are copies of the title 
line of each entry. The current format of the realization 
volume allows a maximum of 9999 lines. The index is not 
considered when determining the number of lines contained in 
eae VOlume. 

Each line of the volume must conform to a specific 
format, of which there are ten: Primitive Title line, 
Comment line, Calc line, Attr line, Call line, Include line, 
If line, Begin Text line, End Text line, and Text line. The 
title line begins with an S or H, to denote hardware or 
software, and contains the name of the primitive. As used in 
the current version of the design system, it is the most 


important of the lines found in the volume, providing the 
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correct format for generating each of the entries contained 

in the intermediate problem description. The calling 
arguments, selection criteria, and attributes of the 

primitive are enclosed in parentheses following ae name. 

The attributes, which vary depending on the nature of the 
primitive, define such parameters as power consumption, 
latency, and chips used. Any or all of the arguments, 
Bemeection criteria, and attributes may be omitted, but the 
commas that separate them must appear. The comment line is 
denoted by the appearance of the letters C-O-M as the first 
Characters on the line. The text that follows is ignored by 
the system. The Calc line allows the use of global variables 
Pecan the system. An example of its application is the 
counter variable used to total the number of chips used ina 
Memercular design. In the current realization library, this 
variable is called ICN and is incremented by any of the 
hardware primitives that contain integrated circuits. The 
Meer line is Similar to the Cale line, but is used to 

Compute the value of an attribute of the primitive realiza- 
tion. Incl and Call lines are used to invoke other primitives 
meom Within a primitive. The difference between the two is in 
the placement of the output that each generates. The output 
from a Call is inserted immediately following any previously 
generated output. Output from an Incl statement will be added 
MOoetnat of the primitive lists after all other output from 


miesineluding primitive has been produced. The If line 
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mm@emrces more flexibility in library construction by 
allowing branch instructions within the realization volume. 
The Begin and End Text lines are used to reproduce descrip- 
tive lines in the files generated for system output by the 
Formatter module. These lines are in the final category of 
exe lines. The most common use of text lines is for 
assembly code associated with a software primitive. Figure 
5 15 an example of a short realization volume. The 
construction of the volume is further demonstrated and 
clarified in the text of the analog to digital converter 


realization developed in chapter four. 
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Weeecs.S 
Ve2e23COM 
V2224C0OM 
V2e225COM 
V2226COM 
V2227COM 
V2228C0M 
V2229C0M 
V2230C0M 
V2231COM 
V2232C0M 
V2233C0M 
V2234C0M 
mee55C0M 
V2236COM 
wee 57 COM 
V2238C0M 
V2239ATT 
V2240CAL 
V2241 INC 
V224eCaL 
V2243C0M 
V2244COM 
V2245COM 
V2246COM 
V2247C0OM 
V2248C0M 
V2249TF 
Vee50INC 
V2251C0M 
V2252COM 
V22S3BEG 
V2254 
wee] 
V2256 
mec] 
V22SBEND 


AMPLE Coleco 0, 52 10,10,—-17,197,19,eeeer 2e58) 
THIS IS A COMMENT DESCRIBING S.SAMPLE. Pi AND 
P2 ARE THE DUMMY ARGUMENTS OF S.SAMPLE. 0 AND 8 


ARE THE MINIMUM AND MAXIMUM VALUES OF THE ACTUAL 
ARGUMENT REPRESENTED BY P1, 0 AND 5 ARE THE 
CORRESPONDING MINIMUM AND MAXIMUM FOR THE ACTUAL 
ARGUMENT REPRESENTED BY Pe. THE 10,10 INDICATES 
ieee oT ORAGERAND IchMe AIMRIBUTES FOR THIS ENTRY 
ARE EACH OF VALUE 10. THE #17 INDICATES THAT THE 
VALUE OF THE EXTERNAL ATTRIBUTE IS CALCULATED ON 
LINE 2259 (2222-(71/7)=2239). THE 18 INDICATES THAT 
LINE 2240 (2e22e+18=2240) CALLS FOR CALCULATION OF 
SOME SGLOBAL VALUE. JHE 19 INDICATES THAT LINE 2281 
GieeormoR The [NGCRUISTON OF SOME OTHER REALIZATION, 
2222 AND 2258 ARE THE FIRST AND LAST LINE NUMBERS 
Ce eS ea 1 ZAeON . 


R EXTERNAL DUMMY 1 7 
Se LOlTEVT RODE Va + 1 
Grd. SL o@onMPLE (0,0) 


Peo eORMPIOKEE (CTOIEVT) 


* 


Heme oOSAMPILE AND S2SAMPIRREE ARE ~TwO GTHER ENTRIES. 
ALSOSAMPLE HAS TWO DUMMY ARGUMENTS, BOTH ASSIGNED 


eevee 2eROn IN THis GASES SeSAMPIAREE HAS ONE DUMMY, 
cep eauAt eG THe GEOBAL TOTEVI FOR THIS CASE. 
DUMMY 1 26GB. 2 SKIP e 
me HeoamerOUR ( } 


INCLUDE H.SAMPFOUR ONLY IF DUMMY1t [IS LESS THAN TWO. 
BeGiNe tee PEXT BLOCK AFTER THE NEXT LINE. 

oN STEXT 

EVES ha tNGeliN THiS SLOGK BETWEEN BEGIN AND END I[S 
CORE tOetme OUTPUT EISTING EXCEPT FOR DUMMY 
ARGUMENTS AND GLOBALS WHICH ARE REPLACED BY THEIR 
ACTUAL ARGUMENTS OR VALUES. 

EX T 


Figure 5. Sample Realization Library Entry 


[Ref. 8] 
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Pie oo eer NGRMETAODS 


With the advent of the microprocessor, a new and 
interesting device for the implementation of digital filters 
Mas been created. Increases in word length and computational 
Seeedewill encourage more complex filtering applications 
Using this device. It is therefore important to investigate 
M@eerecasibility of producing digital filter realizations 
using the automated design system. The fundamental require- 
Memes that the application places on the system can be 
determined using the Intel 8080 microprocessor library 
@mmrently available. Despite the filtering limitations 
presented by its eight bit format, the deficiencies present 
in the current system can be identified and corrected. The 
meeulting improvements will provide a more useful design 
tool and make possible the generation of satisfactory 
solutions to more demanding problems when additional 
realization volumes constructed around advanced 


microprocessors are added to the library. 


fe LRANSFER FUNCTION 
Miemescemeral form of the digital filter transfer function 


is al =) -n 
- ap + ay Z + a5 Z tometer a Z. 
CZ) ==). °° ©. .-2 
ite Ze +t b. Z + .., + b Zz 
il; Z n 


oa 





where ay and b, are real coefficients and n is the maximum of 
the orders of the numerator and denominator polynomials. 
Figure 6 is a flow diagram representation of this transfer 
Mmimetion., There are a number of similar structures that can 
be used to describe the basic filtering operation. Those in 


which the real coefficients ay and b, are multipliers are 


a0 
referred to as direct structures. Nagle and Nelson [Ref. 10] 
describe four direct structures and the equations associated 
with them. 

ie rarst Direct Structure 


Mine first direct structure is derived from the 


equation 


mae coefficient DS Pm qudlisrouone.s Derinane Che imput to the 
filter as X@z) and the output as Y(z), the previous equation 


Can be rewritten as 


Me Ce . etyyt es 6b. 27h) 
baz, } . i oe at : 
1=0 = 0 


This can be redefined using an intermediate variable M(z), 


meen that 
Y(z) Mz) .,™ mee, -i 
Meee i Ee 
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‘The following relationships can now be derived: 


and 


mGz) iL 


er 


Rearranging these equations provides expressions for the 


maput X(zZ) and the output Y(z) in terms of M(z): 
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Because the input is a known quantity, it is desirable to 
Semrye an equation for M(z) in terms of X(z). Manipulation 


of the previous equation for X(z) gives 


n : 
Piez = X(z) - ¢£ dD; z M(z) 


Because the calculations necessary to produce the filtering 
function are to be performed in real time, the inverse 
z-transform of M(z) and Y(z) is taken to provide a time 


mem Solution of the first direct structure. The results 


are: 
n 
Mier = x(k) =— £ bz. m(k-i) 
cet ie 
and 
n 
vik) = £ as m(k-i). 


The flow diagram of the first direct structure, of order 
two, is shown in Figure 7. 
fmeeeccond Direct Structure 
The transpose of a digital filter structure is 
mermed by reyersing the signal flow in ali of the branches 
of the block diagram. Its transfer function is the same as 


That of the structure from which it was derived. The second 
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Perect Beeuerture 1s génerated by applying this principle to 
the first direct structure. The corresponding time domain 
solution also implements the same filtering function, but an 
additional difference equation is required. The equations 


for the second direct structure are: 


p, (k) Ps 4, k-1) it a. x(k) - Ds Cy (k)); He = eve ois. iy tle 


pbk) a_x(k) - Diy (k) 


y(k) = apx(k) + p,(k-1). 


The second direct structure, of order two, is illustrated in 
the flow diagram of Figure 8. 
fee lLhnird Direct Structure 


To obtain the third direct structure, the expression 


memrewritten as 
hn _s n a: 
meee CU. zd) C= XC z)C & ha. zt), 
iz=o * izo + 


The output expression is then 


n 3 n <4 
mG@z) = a, 2 CZ) SY. D. 72 OA 
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Once again taking the inverse z-transform, the time domain 


solution is derived: 


n n 
Pom) = ff a. ita des 8 


Figure 9 illustrates the flow diagram for the third direct 
Mmemeweture of order two. 
HM Fourth Direct Structure 
The fourth direct structure is derived by taking the 
transpose of the third direct structure. The time domain 


representation is: 


Py 6k) = 0) ies r, (k-1) 
q, 6x) = ar, (k) 
More? = =-b r.(k) 
n reyr 10, 
q. (Kk) = a.0Q (Kk) + 434, 6k-1), ht eae g le 


rv. (k) = -bry(k) + r.4,(k-1). 


Spee tlow diagram for the fourth direct structure, of order 


[m7o, 1S) shown in Figure 10. 


Pe oECOND ORDER SECTIONS 
Beene eft the structures described applies to the general 
Pease in which the filter transfer function is of order n. 


The filtering function can be successfully realized using 
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eyo. these methods, provided that n does not exceed two. 
For transfer functions of higher order, coefficient 
sensitivity becomes a significant factor in accurately 
generating the filter output. This phenomenon is created by 
the effects of quantization errors in the representation of 
the transfer function coefficients. The filter structures 
and their related difference equations assume finite 
precision in the filter parameters. However, because the 
word length used by microprocessors is fixed, the precision 
with which the filters can be realized is limited. This 
means that the coefficient representations are not exact, 
resulting in a set of poles and zeros for the realization 
that differ from those desired. The frequency response of 
the actual filter will therefore be different from the 
design specification. In the most extreme case, the filter 
will become unstable if one or more of the poles are moved 
outside of the unit circle [Ref. 11]. As shown by Rader and 
Gold [Ref. 12], high order filters can be redefined as 
combinations of first and second order sections in order to 


memimize these effects, 


fee COMBINATIONAL STRUCTURES 
1. Cascade Structure 
Bees first Combimat1onal structure to be considered 
1s the cascade form of the digital filter. It is illustrated 


in Figure 1l(a) and is obtained by factoring the transfer 
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Sterion into a product of second order sections. This is 


represented by [Ref. 13] 


M 
meze= Tf H.(z) 
isl. 7? 
where 
-1 -2 
ENO al ea a 
H. (2) = = = 
Aly ake Di; Z + Dy. Z 


ad 
wey 921) 12 sen aly 21 


The first element in the cascade receives its input from the 
Signal to be filtered, while each of the remaining elements 
Bets upon the output of its predecessor. 

In z-transform notation, the cascade structure output 


me written 


LaGz7)) -= Hy 6 2) Yay (z) 


il 
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apeGasGaqge Structure 


(z) 





Dm coarallel Structure 


meeure L115 Biock Diagram of Fulter Structures 
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where 


Y, 62) H, (z)X(z) 


ae Z) hele@z Namen ( 2) for 1=2535...,M-L. 

al al i-1 
The equations necessary for a real time implementation are 
formed by taking the inverse z-transform of these expressions. 
For a system of second order modules, the general equation 


memmeeach section is 


fon) = AaguYy_16™ a Ai MYy_y62-}? 3 An uYyiy6n-2? 
~ boyy n-1) ~ boyy (n-2) 
where 
y, (n) 2 ay x(n) a a, 7x (n-l) = Ao1x(n-2) 
- byyyz(n-L) - bay, (n-2) 


y,(n) = 4,593_,6™ + 455¥3_,6a-)) < a5s¥s_1(n-2) 


- by ,y,(n-1) = boy, (n-2). 
The subscript i corresponds to the number of modules needed 
to realize the desired filtering function, and is equal to 
In/2], where n is the order of the denominator polynomial. 


For the case where n is odd, one of the modules will be of 


mraer one. 
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2 Parallel Structure 
The parallel form of the digital filter is the 
second structure to be considered, and is illustrated in 
Figure 11(b). It is formed by factoring the transfer 


mem@etion such that [Ref.,. 14] 


M 
meee = 5 HA, C2) 
fie= 
where 
a.. + a.. go 
2 O1 TAL 
moz) = =I = 
act Dos Zs ct DA Zi 


Ape 1S 271 + a Pi Tee) aS + 4 aa 
_ 0 a 2 3 mn 

FZ) - a =) a8 lh 

ih ae DS Z, = dD, Z + Db, Z. + Dy Z, 

Can be redefined as 
=i = -1 -2 
+ 
ye Ol eM 0 02 2 7 820 7 
= -? =a -? 
io SF Dia Zi + Pal Zz j. ss Di5 Z <r Do Z, 


Fach second order section receives the same input. The 
Sutput of the system is found by summing the outputs of each 
of the parallel elements in the realization, and is 


represented by 
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facez) = 2: 


ie 


Hee z) Xx C2)... 
; 1 


The corresponding time domain equation is 


y(n) = y,(n)ty,(n)+...ty, Cn) 


where 
y,«n) = Ay x(n) 7 a, x(n-l) - biyy(n-1) - Do y(n-2) 


yn) = Ay yX(n) 1 a, 9x(n-1) - Do yfn-1) - Daoy(n=2) 


: + ‘ - - : + > : - A 
y;(n) ay,x(n) a,.*(n I Deiytn Aly) Ds oy(n DY, 

The equations presented for the parallel and cascade 
structures are sufficient to develop the high level problem 
description necessary to realize microprocessor-based 
digital filters using the design system described in the 


previous chapter. 


Meee DESIGN CONSIDERATIONS 

The current implementation of the design system 
generates the control code governing execution of the 
MemliZation using a polled monitor strategy. Pollack 
LRef. 15] argued that this approach was not adequate for 
general controller applications and proposed the introduction 


of an interrupt driven monitor as a user selectable 
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Mieermative. The existing realization volume has the ability 
to recognize an interrupt, but its immediate response is 
iemaced to the setting of a corresponding condition flag. 
The status of the flag 1s examined and the appropriate task 
executed as required on a time schedule basis which conforms 
fo, the structure of the other tasks. The theoretical 
development necessary for implementing an interrupt driven 
monitor has been completed [Ref. 16] and will make possible 
the realization of a wider variety of control devices. 

The polled monitor is the preferred strategy for digital 
filter implementations. Figure 12 is a flowchart represen- 
Tation of three operations basic to the filtering function. 
Fach of these tasks is performed once for each sample taken. 
The interval of time which elapses between successive 
executions of each contingency/task pair is determined by 
the sampling interval required for the signal being pro- 
Meeeead., Only one execution is allowed during that time. 
Therefore, the generation of the filter function is periodic, 
with one output calculated each period of the sampling 
frequency. Because the polled monitor results in the 
execution of the contingency/task pairs at regularly 
determined intervals, it is well suited for applications in 
Mereh periodicity is required. 

Peeesoampling Frequency 

The Nyquist criterion determines the rate at which 


the signal to be processed must be sampled in order to 
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accurately produce the corresponding filter output. The 
amount of time available to perform the required computa- 
tions is then equal to the elapsed time between successive 
Seameres. Ihis, in turn, is equal to the period of the 
sampling frequency. This means that the total time required 
Memeest tor all contingencies and execute their corresponding 
tasks must be less than or equal to the sampling period in 
order to produce a successful single processor realization. 
if this requirement cannot be met, the design system will 
partition the problem into two sets of contingency/task 
Pairs and attempt to generate a solution using two pro- 
cessors. In this case, the contingency/task pairs assigned 
to a given processing element must be able to complete 
Seecution during one period of the sampling frequency. The 
maximum frequency that can be filtered by a given realization 
is therefore limited by the speed with which the filtering 
Seeeeorathm can be executed, which in turn is a function of the 
meeroprocessor family used to generate the realization. 
fee Structure Selection 

The parallel and cascade structures are both readily 
adapted to microprocessor implementation. Because both 
methods are combinations of first and second order sections, 
it is a straightforward procedure to implement each section 
as a contingency/task pair. The natural partitioning that 
exists in each of these structures is well-suited to 


implementation using the automated design system and is 
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meaty Compatible with the partitioning performed for 
multi-processor realizations. 

G@ne selection criterion to be considered when choosing 
a filtering structure concerns the delay time incurred in 
meeawecing an output from a specific input. The delay time 
Will be a constant value for the parallel realization, 
Mee@ependent of the number of second order modules required 
mempertorm the filtering operation. This will correspond to 
an interval equal to twice the period of the sampling 
frequency, assuming a multi-processor realization. The delay 
time associated with the cascade implementation is dependent 
Meemethe number of second and first order modules contained 
Mmeecie realization. The delay time will then be equal to the 
number of modules in the cascade times the period of the 
memes frequency. Therefore, as the order of the filter 
increases, the number of second order modules required to 
implement it will increase, causing a corresponding rise in 
the filter delay time. These relationships are illustrated 
maeraeure 13, 

ieee @emre laeeqd tO ordering Of the second order 
sections within a cascade structure are not of concern for 
Eresapplication presented in this thesis, but in practice, 
Consideration must be given to the changes that can occur in 
The limit cycle and quantization noise properties of the 
filter as the second order modules are rearranged in the 


Cascade. [Ref. 17] 
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Applications problems are important to the development 
of the design system because each one possesses specific 
requirements for implementation. The testing of widely 
varying realizations will indicate those areas in which the 
System is deficient and identify hardware and software 
primitives that are required but abset from the available 
library. The digital filtering problem revealed several 
shortcomings in the current implementation of the design 


System. 


A. PARALLEL PROCESSING 
1. Single Contingency/Multiple Tasks 

The parallel filter structure presented in chapter 
three consists of a finite number of elements that act 
Semeurrentliy on the same data. The output of each processor 
is then passed to a summing element which generates the 
filter output. In terms of the original problem model, this 
corresponds to a single contingency (the availability of 
data for processing) and multiple related tasks (each of the 
parallel processing elements and the summing junction). The 
design system implementation, as well as the specification 
language CSDL, provide no means for expressing such a 


problem directly. 
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meee FORK Construct 
Representation of the single contingency/multiple 
Beeerstructure subject to the limitations of the design 
Manmeuage can be solved by creating a new construct. This 
statement has been added to the syntax of CSDL and is called 
Mere. FORK allows the concurrent execution of two tasks 
followed by a third, which in turn requires data generated 
by the first two. The structure has the general form 
If MAIN then begin 
MemEORK(TSKL ,TSK2) 3 
Ber 1SK3: 
End 
Memen 1s stated in the syntax of CSDL as 
When MAIN: TIME Do FORK(TSK1,TSK2,TSK3). 
TSK1 and TSK2 are the procedures to be executed in parallel 
and are referred to as the "forked" tasks. Because it 
mempbines the results of TSK1 and TSK2 in a predetermined 
meeiton, TSK3 is termed and "joined" Poe The CSDL 
Megemession of the FORK construct is not sufficient to permit 
Miesuse Of the single contingency/multiple task structure in 
design specifications because no translator exists to 
Beoauce the intermediate form of the instruction. However, 
it is possible to manually generate the required high level 
description of the FORK construct and from it determine the 
intermediate representation. To do so requires that a set of 
"dummy" contingencies corresponding to each of the tasks be 


generated. These contingencies test unique flags that are 
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set by an additional contingency/task pair which controls 
the execution of the composite FORK structure. Figure 14 is 
a general CSDL representation of the FORK construct. 

3. Timing Requirements 

The FORK construct makes possible the implementation 
of design specifications incorporating single contingency/ 
multiple task algorithms. The determination of the timing 
constraints associated with the "dummy" contingencies it 
generates presents a new problem. There are two potential 
solutions for consideration. Each of these can be related 
to the single and multiple processor realizations that the 
design system is capable of generating. The former is the 
more straightforward derivation and will be presented first. 
The contingency and task names of Figure 14 are used to 
Meovide Clarification. 

Implicit in both developments is the assumption that 
the user will only specify the time constraint associated 
with the test for the contingency MAIN. The timing require- 
Semes Lor the dummy contingencies Cl, C2, and C3 must be 
determined from this specification. The scheduling of the 
test for MAIN determines how often the concurrent processes 
must be tested. Therefore, the scheduling interval of this 
Contingency will also be used for Cl and C2. This value will 
De denoted as p. The scheduling requirement for contingency 
C3 is related to the availability of data from tasks TSK1l 


eaa §SK2. The worst case condition, that is, the minimum 
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Contingency List; 
When MAINst1 SE Do FORK 
When Cl tte fe Oo TSK! 
When C2 :te SE Do TSKe 
When C3 3:t3 SE Do TSK3 


Procedures, 
Function MAIN 
Binary MAIN,1; 
If GO=t Then MAIN: =1;3 
Exit MAIN 


Function Cl 
Binary Cia; 
If VAR1I=1 Then Cis=1; 
VAR13:=0; 

exo =Cl 


Function C2 
Binary Cauet, 
If VAR2=0 then Cez=l; 
VAR23=0; 

ex3t Ce 


mumetion C3 
Binary C3,1; 
MieOONet=1l [,and. DONEeC=!1 then C33=1-; 
DONE1:=0; 
DONEC 3 =0; 
exit C3 


PEe@emlinwmre obi Ltstame of FORK Construct 
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Task 


Exit 


Task 


Exit 


Task 


Exit 
Task 


Exit 


TSK1 

(generate output, OUT!) 
DONE1L:=1; 

TSK 1 


TSK2 

(generate output,OUT2) 
DONE2: =e; 

TSKe 


TSK3 


OUT3:=OUT1+0UT2; 
TSK3 


Figure 14. (continued) 
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fieerval at which C3 must be tested, exists when Cl and C2 
are always true. If these contingencies are scheduled for 
testing every p seconds, new data will be generated by TSKl 
and TSK2 with the same frequency. This implies that TSK3 
must be executed at intervals of op seconds as well in order 
to maintain proper sequencing of the device. Therefore, the 
meee COMStraint specified for contingency MAIN will also 
Bemeeplicable to Cl, C2, and C3. An important fact that must 
be considered in reaching this conclusion is that contin- 
memeres Cl and C2 can only occur as often as contingency 
wee similarly, C3 can only be true when Cl and C2 have 
been true as well. The timing diagram associated with the 
scheduling of the FORK construct for the single processor 
BeeerZzacion 1S illustrated in Figure l5. 

The scheduling for the multiprocessor realization is 
not as readily determined. As in the serial case, contingen- 
cies Cl and C2 must be tested at intervals of p seconds, the 
iomeeme specification of the MAIN/FORK contingency/task pair. 
Mmemseheduling of C3 is made difficult by the possibility 
that each of the contingency/task pairs may reside in 
different processors and no synchronization exists between 
these elements. If C3 is tested with the same frequency as 
Meena C25 it is possible to "lose" a set of data from TSK1l 
and TSK2. This occurs when these tasks generate data for 
use by TSK3 near the beginning of their respective executions 


and TSK3 processes the data near the end of its computations. 
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If TSK1 and TSK2 begin another processing cycle and produce 
new data before TSK3 has completed its calculations on the 
previous set of data, its output will be determined using 

the second data set generated and the first values will never 
be processed. The timing diagram for such an occurrence is 
shown in Figure 16. 

There are two possible solutions to prevent the loss 
of a data set. The first is to decrease the interval of 
testing for contingency C3. Checking for the presence of data 
from TSK1 and TSK2 with the correct frequency will ensure that 
it 1s processed by TSK3 almost immediately after it is 
generated. This solution results in a new question to be 
answered: how must of a decrease in the scheduling interval 
is sufficient to prevent data loss? The answer is not 
general in nature because the timing specification required 
will be dependent upon the particular application being 
considered. The incorporation of such a construct in the 
design system would require a more detailed analysis by the 
design engineer. This is contrary to the original goals of 
the system and is therefore unacceptable. An alternative 
approach which preserves application independence is to mark 
the data as it is produced by TSK1 and TSK2 and save it for 
further processing by TSK3. This implies the creation of a 
queue in which the data is placed to await action by TSK3. 
The implementation of this proposal requires two queues, one 


each for TSK1 and TSK2. Each time a task generates new data, 


60 





UOTIELUSWSTdWT JOnNAZSUOD YYOJ dosssooudtjztTny uT SSOTJT PLEPGQ “OT dsuUNnBTY 


GG (ANG ee Aik eaniese) 
= eel 1G) 


EL, C0, C0, ENSL +YSeL 
€9D <AduasduTt zUOD 


¢€ Jossa00dg 
iro ae a ies 


éG PEC TZ bed 
a | 204) Z | 20, ZNSE 4SEI 
eae ee aes co. > AdUSsuTIuE7 


Z JOSSsso00dg 
iia eaten ad Wi peed 


ct P1ed TL P2eq 
La 188. a, 168. TEASE OS SIEZ I 
teceainel ve ae To :Aouesut UO) 


T dOsssa0e0dg 
7 eS Peay ee 


bak 





a counter associated with its queue is incremented to record 
the entry of a new element in the service line. The 
Pemeingency C3 then tests the value of each counter. If both 
pmemerecater than zero, the data that has been in each queue 
the longest is removed for processing. The process is 
illustrated in Figure 17. 

The creation and maintenance of these data queues 
solves the problem of scheduling C3. The contingency need 
Smeyepe Lested at intervals corresponding to p. If the queue 
Peumcers are greater than zero, data has been produced and 
is ready for processing. The criteria upon which the 
existence of a true condition for contingency C3 is based 
Meemeneretore been changed from the generation of output by 
MS4l and TSK2 to the presence of data for processing in 
Bieir respective queues. This approach to scheduling of 
the FORK construct does not necessitate the maintenance of 
separate algorithms for the serial and parallel case. The 
Meenot data queues to insure accurate output from a device 
is Suitable for single processor realizations as well. 
Because the execution of contingency/task pairs is synchronized 
fem such designs, data is processed by TSK3 following its 
Beneration by TSK1l and TSK2. The output of each of these 
tasks is stored in its respective queue and is then removed 
mem processing by TSK3. Therefore, the number of individual 
entries in either queue will not exceed two. The CSDL 
specification of the modified FORD construct is shown in 


Figure 18. 
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Contingency List, 
When MAINstl SE Do FORK 
When Cl mie oe Po TSK! 
When C2 ite SE Do TS8Ke 
When C3 3:t3 SE Do T[SK3 


Procedures; 
Contingency MAIN 
Binary MAIN,1; 
If GO=1 Then MAIN: =1; 
Exit MAIN 


Contingency Cl 
Binary Ci,1;3 
Ptmvekt=teethen Cic=1; 
VAR1:=0; 

exyt Ci 


Contingency Ce 
Binary Ce,1; 
Tf VAR2=0 then Ce3=i; 
VAR23:=0; 

exit Ce 


Continaency C3 
Binary Ghea7 


Meee ote and. on ke.gt.0 then €33=1; 
exe C3 


mieteewlo. Coll basting, Modified FORK Construct 
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Task FORK 


Exit FORK 


Task TSK1l 
(generate output, OUT!) 
Q¢12)3:=OUT1; 
CTRiIs=CTR1I417 

Exit TSK1 


task TSKe 
(generate outout,OUT2) 
Q(22):=0UTe; 
CTR23:=sCTReE+1; 

Exit TSK2 


Task TSK3 
QC11)3=QCie); 
Q(el1)3:=QC22); 
OUT3:2=Q0(11)+@QC21); 
CTRiI:sCTRi-1; 
CTRessCTRen-l1; 

Eixayt FORKS 


Pagume. 18.) eCoontinued ) 
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Eee ULI EPLE REFERENCES TO I/0 VARIABLES 
ie Hardware Binding 

The organization of the Intel 8080 realization 
volume as implemented in the current design system is based 
on the concept of hardware binding developed by Matelan 
[Ref. 18]. It specifies that all of the software required 
for a particular implementation must be generated first, 
followed by the hardware necessary to support it. The 
present version of the design program modifies this concept, 
binding the needed hardware to each individual software 
realization as it is generated. This results in the 
duplication of interface hardware each time a given external 
variable is referenced. When the design program encounters 
an input or output primitive, it first generates the 
necessary software followed by he corresponding support 
hardware. The variable name is not compared to previously 
generated I/O interfaces. Therefore, multiple sensing of an 
input or issuing of an output will cause the generation of 
an input/output port for each reference. 

Me soymbol Table Listing of I/0 Variables 

To prevent redundant generation of interface hardware, 
a record of each I/O variable must be kept which stores each 
Variable as it is defined by the user during the intital 
design specification. This can be done very effectively by 
establishing a symbol table which will maintain the status of 


each input/output variable and the hardware generated for it. 
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Rather than producing the hardware required from the software 
primitives representing each I/O operation, it will be 
memeeaced aiter the user defines the external variables in 
the environment section of the design specification. This 
feature will be incorporated in the input module being 


@eveloped for the system. 


C. ANALOG TO DIGITAL INTERFACING 
1. Processor Controlled Conversions 

The realization library available in the current 
implementation of the design system allows communications 
between dedicated microprocessor systems and external analog 
Signals using a Burr-Brown ADC82AG eight bit analog to 
digital converter. Synchronization among contingency/task 
pairs is a critical factor in accurately generating digital 
filter output. Therefore, the sampled data input must be 
available from the a to d converter on a periodic basis as 
requested by the main processor. The converter available in 
the Intel 8080 realization volume is controlled by an 
independent clock and is therefore not suitable for such an 
@eplacation. The addition of an analog to digital converter 
Beatable for digital filtering applications is therefore 
warranted. 

The device chosen for inclusion in the realization 
library is the Analog Devices AD570 and is an appropriate 


selection for several reasons. The most significant of these 
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is its ability to perform a conversion when directed by an 
outside signal. Connections are provided for both a data 
convert signal and an associated output to indicate the 
Semeletion of the conversion operation. The circuits 
required for interfacing the device to the Intel 8080 data 
bus are provided in the manufacturer's specification sheets 
(Ref. 19]. The device is also capable of accepting unipolar 
SeeoapoOlar data input. 
fe Gibrary Entry Development 

The generation of the library entry for this device 
is a simple task. Because the use of hardware in any 
realization is controlled through software primitives, a new 
software entry was generated as well. This new primitive is 
called S.ANAIN. User specified parameters for the primitive 
are determined by those required for inclusion of the 
necessary hardware realizations. These consist of the name 
Of the output signal, the maximum and minimum values that it 
can assume, and the number of bits in the digital represen- 
tation. These parameters represent the data required for 
correct generation of the a to d converter hardware. The 
3dB bandwidth of the input signal is also required in order 
to produce a buffer amplifier to the input of the device. 
The high voltage specification is used to determine the gain 
provided by the buffer amplifier. The interface hardware 


listed in the specification sheets is not required for 
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[miemementation by the design system. The inclusion of the 
primitive S.SENSECOND provides a suitable interface for the 
Sevice, 

The a to d converter is listed in the hardware primi- 
tive H.ADC2 and contains the pin connections necessary to 
implement the device. The low voltage limit is evaluated to 
Seer mine if the interface for bipolar input is needed. The 
Giebe line for both of these primitives is copied to the 
index of available realizations that is at the beginning of 
the volume. The attributes of time, storage requirements, 
E@oaeinciusions are listed as well. A complete listing of 


Seemeprimitive 1s contained in Figure 19 and 20, respectively. 


eye 





Seanain Msn annvl, 01070, 20,-e5, 1,100:,4,7, 5916, 3847) 
com primitive to define orocessor controlled analog input 
pomeryst = wmoutrhi volt limitrelo volt limit,s3db rolloff: 
com bitsrs,voltage limits,bw limits: 
com time,stor,ext,calc,incl,addrs 
com added by mere. heilstedt, may 1983 
name $na001-5na030 
begin stext ; 
ranalog inout channel for signal <sig>, 
srPange <h> to <I> volts 
moodbp rolloff at <b> khz 
endtext 
incl h.conn-al (<$na006>,<$na007>,QGgnd,<sig>::) 
com select again for buffer amo to match inout range 
if <I> .It. 9 skip 4 
yr <h> le. 10 skio 6 
if <h> .le. 25 skio 8 
mtiescn> ete. 50 skio 10 
com set buffer amo if bioolar output 
mtesl> .ge.e “5 skio 3 
meee > .9e. “ie skio 5 
feel? .@Q. -25 skio 7 ; 
com gain 1.0 (written 10) for input range 0,+10;7;+-5 volts 
incl hebufframp (<$na006>,<$na007>,<5na005>,10,<b>3:) 
skip 5 ; 
com gain 2.5 (written 25) for inout range 0,+25;+"12.5 volts 
incl hebufframo (<$na006>,<$na007>,<5naN05>,25,<b>:3 ) 
skio e : 
fom dain 5-9 (written 50) for inout range 0,+507+=-25 volts 
incl hebufframo  (<$na006>,<$nal007>,<3naN005>,50,<b>::) 


incl h.adce (<$na005>,<sigr>,<h>,<1>:8:3) 
call sesensecond (<sig>:8,128) 
com 


Figure 19. Analog Input Primitive 
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neaoce (anrouteis | f0rO. per er SO4G , 5899) 
com erimitive to define 8 bit orocessor controlled adc 
com list = inputroutputshi volt limite lo volt limit: 
com bitsslatrowrrchiopsrcalcrinel,adars 
com added by me re heilstedt, may 1983 
name $na031-5na060 
incl h.trimpot (gnd,<$na041>,<in>,2002) 
begin htext 
a/d converter, 8 bit 
device is analog devices ad5/0- ic <icn> 
connections: 
pin 1 = Miae< 


DIN ee = <out>(0) (1sb) 

pin 3 = <out>(1) (2lsb) 

pin 4 = <out>(2) (31 sb) 

0oin 5S = <out>(3) (41 sb) 

DIN 6 = <out>(4) (53 sb) 

pin 7 = <out>(5) (61sb) 

pin 8 = <out>(6) (71sb) 

pin 9 = <out>(7) (msb) 

pin 10 = +5 volts 

oin il = conv (blank and .not. convert) 
pin le = -15 volts 

pin 13 = <$na041> (Canalog input) 
oin 14 = gnd (analog common) 


endtext 

men ii> it. 0 skip 6 

com if unipolar inoutsy make direct connections 
begin htext 


pin 15 = gnd (bipolar offset ) 
pin 16 = gnd (diqital common) 
endtext 
skip 11 
hegin htext 
pin 15 = <$na049> (bipolar offset) 
pin 16 = <$naN48B> (digital common) 
endtext 
inc!) henand4 (O,1,+5v-<tnaNs4B>, <G$naN04S?:) 
incl] h.einvert (<$na045>,<$na046>:: ) 


inc) hediodewsw (<$na046>,<$nal4/7>2: 
incl hediodewsw (<S$na047>,<$nal48>?: 
incl hediodewsw (<S$na048>,<$naN49>22 
incl] h.eresmfaqtrwt (-15v,<5na048>, 50000 
begin htext 
Sin i? 
pin 18 


) 
) 
) 
ae 


dr (data ready) 
ec. 


endtext 
cale icn=icnel 


Figure 20. Analog to Digital Converter Primitive 
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VY. DESIGN EXAMPLES 


The ability of the design system to produce an acceptable 
realization is demonstrated using a realistic problem taken 
from the article by Nagle and Nelson [Ref. 20]. The 
implementation considered is a fourth order digital rate 


mm@irer with the transfer function 


a S 4 


#+1.77044z 
im 


#0. 3443332 


$3.5368227--1. 8886727 +0. 39750627 


1.04+0.3902u42 7+], 2u2u727 
1.0-3.02828z7~ 


BZ) = 
This expression is redefined as a combination of second 
order modules from which the difference equations necessary 
to generate the filtering function in real time are derived. 
The choice of implementation algorithm for these modules is 
a difficult one because none of the four methods presented 
in chapter three is clearly superior. Each of the direct 
forms consists of the same number of multiplication and 
addition/subtraction operations. However, the direct form 
one algorithm possesses the minimum storage requirements of 
the four candidates and will therefore be used in the 


development of the example designs. 


fe 





Pee CASCADE REALIZATION 
The cascade implementation of the proposed filter 
consists of n/2 or two second order sections. The equivalent 


Expression is 


Meee Hi, (z) : H(z) 
_ «A+. 79062774+1.287227° | 1-2.18827~+1.38322 ~ 
14#1.82327-+0.89592- 14+1.205427++0.4u3827° 


im CSsdL Description 

Despite the fact that the high level design language 
Meee will not be implemented in the input specification 
module, it remains a useful tool in the current limited 
implementation of the system, providing a means of trans- 
lating the problem statement into its intermediate form. 
Therefore, the first step in the development of the desired 
mebter implementation is to express the problem in the 
Seaeax of CSDL. 

Micmraecitarication section on the listing contains 
all of the administrative data associated with the problem, 
such as the designer's name, date of creation, and project 
pele, ft is as follows: 

IDENTIFICATION; 
Deememom: Ma Rk. Heilstedt 


Date: 4Y4-7-83 
Project: Sample Cascade Realization 
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The environment section is more easily determined 
meeerm che completion of the other subsections in the CSDL 
listing. It will therefore be developed last. 

Miew2emeracion of the contingency/task pairs is 
accomplished with the aid of the flowchart shown in Figure 
21. Each of the decision diamonds represents a contingency 
to be tested, while the process boxes associated with them 
are their corresponding tasks. The contingency/task pairs 
and their order of specification in the contingency list are 
therefore determined by the sequence of execution indicated 
mothe flowchart. 

The first contingency/task apir determines when a 
sample of the input signal must be taken. Because this must 
be accomplished in accordance with the sampling theorem, the 
"dummy" contingency EVERY is used to force the generation of 
a sample once each period of the sampling frequency. 
Arbitrarily choosing a sampoing interval of sixty milliseconds, 


me first pair is 


Evemy OO0MS do SAMPLE. 


The second contingency/task pair tests for the 
Peartability of data from the analog to digital converter. 
Because the Analog Devices AD570 takes approximately twenty- 
five microseconds to perform a conversion, a test for the 
Beesence of data at its outputs is required. this is 


accomplished by sensing the value of the single bit input, 
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data ready. The variable DR is chosen to represent this 
Control line. The value of this input is normally zero, but 
will change to logical one while the conversion process 

takes place. A return to zero indicates that the operation 
has been completed. A sample of the input signal is required 
at intervals determined by the period of the sampling 


frequency. The second contingency/task pair is therefore 


When READY:60MS do READ. 


The third contingency/task pair tests to insure that 
the reading of data has taken place and when true, executes 
mem tirst filtering function of the cascade. The time 
constraint associated with this pair is once again equal to 
the period of the sampling frequency because data for 
processing is produced with that interval. The third 


contingency/task pair is then 


When NEWDATA:650MS do FILI. 


The final contingency/task pair in the list tests for 
the generation of data by the first element in the cascade. 
When true, the second filtering function is executed. 

Because the contingency can be satisfied only as often as 
data is produced by the first filtering function, its time 
Sonstraint is also equal to the sampling interval. The final 


contingency/task pair is 
When ONEOUT:60MS do FIL2. 
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imies procedures section of the CSDL listing requires 
Peeopecification of the high level code required to test the 
contingencies and execute the tasks. The former shall be 
developed first. 
Mge contingeney for EVERY is a "dummy" specification, 
used to force the execution of a specified task at user 
Meeimed intervals. As such, it has no high level listing. 
The contingency READY senses the value of the input 
Hine DR. When DR is equal to zero, READY is set to one. 
Because the analog to digital converter will require 
twenty-five microseconds to complete the conversion operation, 
erred Walt corresponding to this period is inserted in the 
contingency to guarantee that sufficient time for the 
Senversion to take place prior to testing for its completion. 
The contingency is 
Contingency READY 
Binary READY,1; 
Waa 2ouS: 
Sense(DR); 
Tt DR=0 then READY: =1; 
Ee READY 
Wie tied Genzringency, NEWDAITA, tests for the value 
amor iae, Called REDDAT, which is set at the conclusion of 
the task READ. A value of one indicates that data has been 
read and is awaiting processing. The contingency is 
Contingency NEWDATA 
Binary NEWDATA,1; 
If REDDAT=1 then NEWDATA:=1; 
REDDAT:=0; 

Exit NEWDATA 
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Gime tinal contingency, ONEOUT, tests for the 
mei lability of data from the first filtering function. When 
Ieee value of the flag DONE] is set to one at the end of task 
Meise data 1s ready for processing by task FIL2. The final 
Contingency is 
Contingency ONEOUT 
Binary oneout,1; 
Beene then ONEOUT:=1; 
DONE1:=0; 
eet ONEOUT 
In addition to generating the filter output, the tasks 
also perform internal operations, such as the setting of 
eee which are invisible to the user. The first task to 
memomplemented is SAMPLE, which is used to produce the signal 
required to initiate a conversion by the analog to digital 
converter. This is done by setting the control signal, CONV, 
memzero, issuing it as an output, resetting its value to 
Mogical one, and issuing it again. The listing for the task 


ies 


Task SAMPLE 


GON .=0 5 
iesue (CONV); 
COny;s=1; 


Issue(CONV); 
Peele SAMPLE 


Task READ senses the value of the input variable as 


determined by the analog to digital converter. It must also 
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set the contingency READY to zero and the data available 
fetes RKEDDAT, to one. The high level description of READ is 
Task READ 
READY: =0; 
Sense(X); 
REDDAT:=1; 
Exit READ 
The tasks that generate each of the filtering 
operations are similar in design and will therefore be 
meweloped together. FILL resets contingency NEWDATA to the 
false condition. The filtering algorithm is executed, with 
the result being stored in a global variable. The flag 
Pomel 1s set to one to indicate the completion of the 
Computation, and the intermediate variables used are updated 
for use in the next computation. FIL2 first sets contingency 
Seoul to zero. It performs the same general calculations 
[eeepe ditrterent coefficients, and issues the result of its 
Gomputations as the system output variable Y. The task 
Meocings for FIL1 and FIbL2 are 
Task FIL 1l 
NEWDATA:=035 
amiga = ( 182.2852 ) = G08 959"M13) ; 
cae ee 7 Wo leg ee oe Mies 
DONE. +L: 
Migs Mad : 


M12:=M1l1; 
Eee FLL 


fee 





Task FIL2 
OME Unk es SNe 
ee a 2054 M 22 = CO S44 38"M23): 
ie eee M22) + (238s 2*M23): 
Tssue(Y); 
Ma oe M225 
M2 2eraiZ 1: 
Beat FIL2 
hanes Che Completion of the contingency and task 
specifications, the environment section can be defined. Each 
of the input and output variables must be listed first. They 
are represented by X and DR, and Y and CONV, respectively. 
The arithmetic variables are the global parameters used 
within the CSDL listing by multiple contingencies and tasks. 
Listed with each variable is a description of the type of 
Signal it is and the number of bits required to represent 
it. The completed environment section is 
ENCIRONMENT ; 
Input: X,8,TTL;DR,1,TTL; 
Meepucs  Y,6,ClILsCONV,1,7IL; 
PwaminmetLle:) Yl, 242y2,24=M11 24 :M13,24; 
ee Meee eee eae M3 62 ae 
M33,24;READY,1;DONE1,1; 
VNEVDAtAS On aOUlT i. 
Rew iA ae 
2. Intermediate Representation 
The generation of the list of primitives and the 
MmepEFL file for the example specification is a line by line 
translation of the CSDL listing. The IADEFL file is derived 


heom the identification section and the contingency list. It 


1s produced by extracting the names of the contingency/task 
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mers and their timing constraints from the former and 
Mestang them individually. In addition to the period of the 
contingency, the user may also specify the maximum time 
allowed for the execution of the Sone enee the maximum 
allowed duration of the task, the global order of the 
contingency/task pair, its priority, and the maximum allowed 
time duration of any timed block within the contingency and 
its maximum allowed duration. The metric upon which the 
design is to be generated is not part of the CSDL specifi- 
Seeron but is included in the IADEFL file. The available 
criteria for design selection are: first successful 
realization produced, realization with least power require- 
ments, and least costly realization. Because it is the only 
one implemented, the first metric was chosen. The design 
identification data follows the listing of the design 
criteria. 

The procedures section of the CSDL listing contains 
the information needed to generate the primitive list 
Specification of the design problem. Each line of the high 
level expression is translated into one or more lines of 
intermediate code, In order to manually generate the 
primitive problem representation, the user must have 
available a listing of the index to the realization volume 
used. Each of the entries in the intermediate specification 
is derived from the title lines of the available primitives. 


The process of translating the CSDL listing into its 
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intermediate form is tedious and repetitious and will not be 
@eteiled here. However, a short example is instructive. The 
line to be translated is a mathematical expression taken 

meeme the first filtering task, FILL. The example statement 


Le 


ae = x—(1.823"M12)-(0.8959*M13) 


The corresponding primitive code is produced by working from 
meme tO left in this statement. Each of the expressions in 
parentheses is derived first. The results of these operations 
are then combined as indicated by the remaining operations to 
attain the desired result. The constants must also be 
specified as variables within the listing. Each variable 

meme used requires an additional statement for the genera- 
mm@eon ot 4a storage location. The primitive list for this 


equation is 


s.fmul (mla,allj,m12:24,24,24) 
Seemu 1. (mlb al 2emils : 24.20.24) 
s.fsub Gianmla.mibe24a24 26) 
s.float Cs alioa mr, dis ano, 5) 

s.fsub Cm, Sicinoma: ct. 2424) 
S.var miluas 2/1) 

S.var Cre 2a 21h) 

S.var Gnd =) 

S.var (m13:24) 

S.var (ma: 24) 

S.var Cinsig: 24) 

s.fcons Calis, OO 0s 0e 24) 
s.fcons Ca 17eO. Ol. 0 Or 21) 


mie TADEFL file, primitive list, CSDL problem statement, and 


design system output are includes in Appendix A. 
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Eee PARALLEL REALIZATION 

The parallel realization of the fourth order digital 
rate filter also consists of two second order modules. The 
Filtering function to be generated is found by taking the 
Meweda! fraction expansion of the original transfer function. 


The resulting expression is 


ie cng eo 


ee eo ie Oe 9 O19 Z, 





¢ 


1 iz 


ewe? sa 56927. 
=) 


ee ee 
1a so a0, Wines 


2 


Memes 1llustrated by the flowchart of Figure 22. 
m weobL Description 
The CSDL description of the parallel realization is 
Similar to its cascade counterpart. Therefore, only those 
areas in which they differ are presented in detail. The 
m@emrification section for the parallel realization is 
IDENTIFICATION; 
Designer: M. R. Heilstedt 
Project: Sample Parallel Filter Problem 
Dace “ieaprit 1983 
The dummy contingency EVERY is once again used to 
fonee the sampling of the input signal according to the 
interval determined by the sampling frequency. The second 


contingency tests for the completion of the conversion and 
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executes the task which reads the digital data input. Each 
of these contingency/task pairs is identical to those of the 
Same name in the cascade problem specification. 

The remaining entries in the contingency list are 
unique to the parallel problem and are used to implement 
m@meerORK construct. The first of these tests for the 
availability of data for processing. A true condition 
executes the task FORK, which sets two variable flags, VAR1 
SemveR2, tO one. Contingencies Cl and C2 test the values 
of these flags. When true, Cl permits the execution of FILI1, 
Eiemtirst Of the parallel filtering functions. Similarly, 
feemsecond filtering slgorithm, FIL2, is executed when the 
contingency C2 is satisfied. FIL1 and FIL2 each place their 
data in individual queues and increment the value of their 
Memeipective queue counters. Contingency C2 tests for the 
existence of each queue and if both are present, causes the 
execution of task PLUS. Once again choosing an arbitrary 
mamerine interval of sixty milliseconds, the contingency 
mest is 

CONTINGENCY LIST; 

EVERY 60 MS do SAMPLE 

When DATA:60MS do READ 
MiiemmbeOr-oUMS ido FORK 
Mine nn Cees UM Sa co) Ei ied! 


When C2 -s0MS dow FE lilZ 
When C3 oO UMSmcdo.fbUs 


eine the terminology introduced in chapter four, FIL1 and 


FIL2 represent the forked tasks and PLUS is the joined task. 
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ime algorithm for contingency DATA is identical to 
feaeesed for its counterpart in the cascade implementation. 
Contingency LOOP tests the value of variable GO which is set 
Memone by the execution of task READ. When the contingency 
is satisfied, it permits task FORK to be executed. The 
Peyeingency is 
Sontingency LOOP 
Bamany. LOOP si: 
hae SS IL te loveyol BLNele! 2 lee 
GO=0; 
Beet LOOP 
Contingencies Cl and C2 perform similar functions, 
using the variables VARI and VAR2, respectively. When true, 
Meepermits execution of FILL and C2 permits execution of 
FIL2. Both of the test variables, VAR1L and VAR2, are set by 
execution of task FORK. The contingencies are 
fentingency Cl 
Baars Ll: 
If VAR1=l1 then Cl=l; 
VARL=0; 
meee C1 
Contingency C2 
Emeilaiyec 2,5 
iar VAR2sweethen C2=15 
VAR2=0; 
esa t. C2 


Monermeeney 3s tests the value of the queue counters 


to determine if data generated by tasks FIL1 and FIL2 is 
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Waiting to be processed. If both queues contain data 
P@erres, task PLUS is executed. The CSDL listing for C3 is 
Cemtingency C3 
Bamooy C3 ,.1% 
Meee O ssand,. CIR2,.gt,0 then C3:=1; 
Eyer C3 
The tasks SAMPLE and READ are identical to the tasks 
already described for the parallel filter implementation 
problem. 
Task FORK sets the flags tested by contingencies Cl 
and C2 to one to represent the true condition. It also 
resets the value of its associated contingency, LOOP, to 


zero. The task listing is 


Task FORK 
E@Or: =O: 
VARL:=1; 
VAR2:21; 

mett FORK 


It is important to note that the contingency/task pairs that 
follow cannot be satisfied if the LOOP/FORK pair is not true 
as well. Because the tasks that perform the filtering 
operations are included in this group, FORK effectively 
Bemerols the execution of the filtering function. 

The CSDL descriptions of tasks FIL1 and FIL2 are 
Similar to those used in the cascade realization. The direct 
form one algorithm is again used to implement the filtering 


mmerctone. Ihe results produced by each of the tasks are 
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placed in individual data queues and their associated 


counters are 


the tasks are 


Task 


Eee t 


eis ic 


pet 


ineremented. The high level descriptions of 


IB SCAB 

Cl y=0e 

aie Cer sii 2)=(04596"ML3); 
Q1B:=(1.264*M11)-(1.748*M12) ; 
Maks = Ml 2 

MU? oak © 

Clue le = C ip abaal:: 

Prod 


Rie 

C2 = 0)s 
M21:=X+(1.205*M22)-(0.444*M23); 
OB iClns / 23M 569-2 2 ) ; 
MZs" =M2 2 ; 

WAZ Saez 1: 

Cine, =e not LS 

|B ae 


Task PLUS adds the data generated by FIL1 and FIL2 to 


Meeamee the output of the device. The queue counters are 


decremented following the output of the filtered signal. The 


Swe description of PLUS is 


Vask PLUS 


Y:21.0-(8.0*Q1A)-(4.0*Q2A); 
Issue(Y); 

Cie =CiRl=1: 

Ci? s=CTRZ=L: 


Bee PLUS 
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Now that the contingencies and tasks for the parallel 
| problem have been defined, the environment section entries 
| can be determined. 

ENVIRONMENT ; 

iepuc: X,8,lTb3DR,1,TTL: 

Pmeput: Y,8,2fTL;CONV,1,TTL; 

Pyemenmette: OQ1LA,24:01B,243;02A,24;02B, 24; 

Ma 2 oe Mia M21 24: 


Zee 7 ia GOe bs VARI 8: 
VARZe LOOP so DATA .8 : 


2. intermediate Representation 
The IADEFL file and the primitive list are generated 
Meemecne CSDL listing as detailed in the development of the 
input specification for the cascade problem. Appendix B 
Gontains these listings, as well as the complete CSDL 
specification of the problem and the output data generated 


Pomeeme desagn program. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


eee «CONCLUSIONS 

mie feasibility of automating the production of 
microprocessor-bDased digital filters has been demonstrated. 
The use of applications problems has been shown to be a 
logical and productive step in the development of the design 
system. The modifications and additions made as a result of 
this research have increased the versatility of the current 
version of the design program and identified areas for 
fmeure study. 

The value of the design system described and tested in 
mms thesis is readily apparent. The time now required to 
manually produce the hardware and software necessary to 
Support dedicated microprocessor systems can be more 
productively spent if the computer can be programmed to 
Beecomplish this for us. Further development of the design 


system is therefore warranted. 


Bee RECOMMENDATIONS 
Ll. Implementation Language 
The current version of the design program is written 
in FORTRAN. This language places rigid requirements on the 
format of the input data files used by the system and 


restricts the organization of the program code. Modularization 
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ene program algorithm into specific levels of operation is 
Mee possible. The use of a language such as Pascal, PL/I, 
or Ada would permit this type of program organization, 
allowing additions and modifications to be made more easily. 
Validation, verification, and testing of the design program 
would also be enhanced. 
2. Validation of Current Program 

The design program installed on the Digital Equipment 
VAX 11/780 was received on magnetic tape from Lawrence 
Livermore Laboratory. Because of a compatibility problem 
between the machine that produced the tape and that which 
read it, sporadic errors occurred in the copying of the tape 
onto the VAX. These errors consisted primarily of incorrect 
Dranch statements, but erroneous variable references were 
also found. These problems were of sufficient severity to 
prevent realizations from being generated when in fact they 
were possible. Therefore, a considerable amount of time 
was spent debugging the design program and in fact became a 
major portion of the effort to produce this thesis. A 
thorough validation of the program remains to be accomplished. 
That which has been done thus far has been performed to 
identify the source of observed errors in program execution. 
Due to the nature of the inconsistencies found it is 
reasonable to assume that those subroutines not checked 


contain errors as well. A line by line manual comparison of 
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Mj-mlascting of the design program installed on the VAX 11/780 
Depend COpy of the code that is known to be correct is the 
Best method for determining the correctness of the current 
implementation. 
fee Realization Library 

The choice of implementation language for the design 
Peegram Will also determine the format of the entries in the 
realization volume. Its current organization makes additions 
and modifications a difficult task. Independent of the 
implementation language and library format is the need for 
an algorithm to allow library updates by the user. This type 
of program would save considerable time in the development of 
the hardware/software database. 

The organization of the library merits further study. 
In keeping with the concept of hardware binding, the user 
specifies his problem in terms of software primitives. The 
support hardware necessary 1s automatically generated. 
Because the design algorithm is intended to search only one 
volume at a time for the needed primitives, duplication of 
hardware between volumes can occur. To prevent this redundancy, 
the creation of a global realization volume of commonly used 
hardware primitives would be helpful. After an unsuccessful 
search of the current microprocessor volume, the program 
would scan the global listing in an attempt to locate the 
needed primitive. Failure to produce a realization would 


Beeur only after an unsuccessful search of both the current 
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and global volumes. The content of the global volume would 
not be limited to modular hardware. Individual components, 
such as resistors and capacitors, could be included as well. 
Miesprinciple of hardware binding is not violated as long as 
user access to these entries is only allowed through software 
primitives. 

The ultimate goal of the design system is to produce 
a physical hardware realization of the problem specification. 
mmouwehn a system, the current program would be the first 
module ina group of three or more that would generate the 
layout for the device program read only memory, with the 
monitor and désign program, and assemble the hardware. The 
program code required to implement such a system would 
necessitate the listing of the software primitives in terms 
of microprocessor operational codes rather than assembly 
language. Therefore, the organization and incorporation of 
op code based volumes in the realization library must be 
mavestigated. 

eeeeecerrupt Driven Monitors 

The theory for the implementation of an interrupt- 
driven monitor has been developed but inclusion in the 
design system remains to be accomplished. The availability 
of such a strategy as a user selectable alternative will 


greatly increase the potential of the design system. 
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5. Applications Problems 
More applications problems are required to determine 
Siesremaining shortcomings in the design system. Signal 
processing implementations involving modulation/demodulation 
techniques, as well as high speed special purpose hardware 
realizations, remain to be tested. The availability of the 
interrupt-driven monitor strategy can be expected to create 
additional test applications. 
oe Documentation 
The documentation currently available is inadequate 
mepermit rapid familiarity with the program. Therefore, 
reference data must be maintained which details the 
development and implementation of the system. This will aid 
subsequent research and provide the basis for a user's 


manual when a commercially useable design system is produced. 


Qu 





AEPENDIX A 


Dole CASCADE REALE ZAT ION 


Siiseappendix contains the input and output files for 


the cascade realization developed in chapter five. 


A. INPUT DATA 
1. CSDL Listing 


Mees 27S the CSDL description of a" 

Bourti order digital 'rate' filter. It is" 

"taken from the article by Nagle and" 

"Nelson which appeared in the February 1981 issue" 
MereCOMPUITER magazine. It is implemented as" 

Mere cascade of two second order sections, which" 
tare expressed using the direct form one algorithm." 


eieN TIFICATION ; 
becmeomer:M. R. Heilstedt 
Datex4-7-83 
Projett: Sample Cascade Filter Probiem 


ENVIRONMENT 3 
ieee ol LL3;DR,1,IIL; 
Ouepmees Y,6,.1L;CONV,i,TIL; 
Pommem@eret eG: Y1,8:¥2,83;M1l1,8;3;M13,383M21,83M22,8; 
Mis ee otha 2 50h M Shon aivi co aie aia: 
DON ae: 


Sey INGENCY LIST; 
mSample signal every 60 milliseconds” 
EVERY 60ms do SAMPLE 
teheeck for conversion completed" 
When READY:60ms do READ 
iweerinst second order filtering" 
When NEWDATA:60ms do FILI 
Tjeneae second filtering" 
When ONEOUT:60ms do FIL2 


his. 





meoCEDURES; 
Contingency READY 
Binary READY ,1; 
Mett 25US; 
Sense(DR); 
iMieeDR=-0 then READY:=1 fi; 
poi, READY 
Contingency NEWDATA 
Binary NEWDATA,1; 
If REDDAT=1 then NEWDATA:=1 
REDDAT: =0; 
econ NEWDATA 


Semermesency ONEOUT 
Eanary ONEOUT ,1; 
fe WONEL=1 then OQNEOQUT:=lL; 
WONE L:=0; 

eet ceONEOUT 


Task SAMPLE 


CONV:=0; 
Issue (CONV); 
Gaiie= |: 


Tsgue(CONV); 
Exit SAMPLE 


Task READ 
READY: =0; 
Sense(X); 
REDDAT:=1; 

Exit READ 

Mase b LT L1 


NEWDATA: =0; 
iets == 16 73M 2)=(0.8959"M13)-; 
reel + C7 906*MI2)+C1, 287 2*M13): 
Soe: =1; 
eset il 2 5 
faeZ > =Mi1:; 

peer LL) 


Macicee L152 
eVEQULs=0: 
Pe s=Yt—( 1. 2054"M22)=-(0.4438"M23): 
aera 2 2 BS M22)+(1.3832°M23) : 
Issuel(Y); 
as 2 =M2 2 ° 
MzZ22=M21: 

perc le ii 
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WeOEFL File 


ie. 


aS oo © © 


2 ® & & & 


oo °o © 


.» ® & & 


= 


=o 2 © 


ss ® & j& & 


= == <= = 


—~ = = = 


W31G0Nd YILTI4 AGVISVI AIdWVS 


SE 24TsE fos ps assy 


clts 


Tits 
peas 


ataues 


NOBUO 
eVepmMaus 
eyep: 
yoeas 

Wa SAS; 


S00 
100 
£00 
200 
100 


0000 0 Dee ee 


sli) 





ae © Ua © Iam © um © ur © Dm © Dam © 2m © Daum © Tuam © Ju © Ju © Juan © Juan © 2ums © 2am © tun © 2uums © un © 2am © Dus © Jum © Jas © uum © Juus © Daum © Jams © Ju © Faas © 2 @ Joa © Ju © Dau © Das © 2m © 2s © un © 2 © Eu © Dan © Da © 2 © 2 © Da © a © 


3. Primitive Listing 


t. generated forssystem KR RR KR KRHA KK KK 
Semain (53535) 

S.start C333) 

te generated forteach KKRAKRKKK KK KKKAKKKKKKAK 
s.every Cegchs :) 

Sevar (each:8) 

t. generated foridata RRR KKK KKEKREKRKEKKKKKKKKEK 
S.O ROC (daeas:) 


s.fixedwait (25:3) 
Ssesensecond (dri1,1) 


$eeQq (dte,der,9c00128,1,1) 
Sejmof (dte,d100138) 
Seasstaqncons(data,1:1,1) 

S$.}o0c (9100132) 

Ss.exitoproc (data,0::) 

$eCONS CJe001T, 921; 1) 

S.var (deel) 

Sevar (data:1) 

S evar, (ite: 2) 

t. generated forinewdata eK K ARK KKKARKK AK KKK KKK 
ScD EOC (newdata::) 

$.eeq . (dt 3,reddat,92c002:8,1,1) 
S.jmof (dt 3,901002:8) 


Seassigncons(reddat,0:1,1) 
Seassigncons(newdata,|!1:1,1) 


seloc (91002:8) 

seexitoproc (newdata,0::) 

$.CONnS (9¢002,031,1) 

Sevar (newdata:!1) 

Sevar (reddat:1) 

Sevar (3¢3:8) 

t. generated fortztoneout RHEEKKKKEKKAKKKKRKKKKKAKE 
S10 0C Comeout > 3) 

$eeQ (at 4,donel,ac0035:8,1,1) 

Ssejmof (9t4,01003:8) 


Seassiaqnconsl(oneout,1:1,1) 
Seassiancons(donel,0:1,1) 


seloc (910033) 

seexitoroc (Coneout,0::) 

S$econs Caevos, 121, 1 ) 

Sevar (oneout:1) 

Sevar (donel:1) 

Svar (dt438) 

te. Qenerated for:samole kaK KK KKK KKK KKK KKKKKK 
Se-Proc (sample::) 


Seassigncons(conv,0:1,1) 
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seosuecond (conv:l1,c) 
Seassigncons(conv,1i:1,1) 
S$-1ssuecond (convil,2) 
Seexitproc (samole,each::) 


Sevar (conv:l) | 

t. generated foriread KRHA KKRKKRRAKAEKRKKKKKKKKR 
Seproc (read: ) 

Seassigncons(data,0:1,1) 

Seanain Cx 107 10,5:8s) 
Seassiaqncons(reddat,1:1,1) 

Seexitoroc (read,data::) 

Sevar Exe) 

t. generated for:fill KRARRERKKRERKAKKKAKK KKK 
S$.Oroc Chi 1s: ) 

Seassigncons(newdata,0:1,1) 

se fmu] (mlbo,be»m132:24,24,24) 

s.fmul (mlarbi,m12e:24,24,24) 

s.fsub (marymiarmlb:24,24,24) 

sefloat C(ix,x38) 

se fsubd C(nlil,ix,mas2e4,c4,2c4) 

s.fmul (yib,a2,~m133:24,24,24) 

s.fmul (viar,al,mi2:24,24,24) 

safedd (yarvylar,ylb:24,24,24) 

s.fadd (yviesmil,yaz:e4,24,24) 


Seassigncons(donel,1:1,1) 
sefassigqn (m13,m123:24,24) 
se-fassign (m12e,m11324,24) 
Seexitoproc (fill,newdata:: ) 


Sevar (mibsed) 

Sevar (m133:24) 

Sevar (mila:2e4) 

Sevar (miered) 

Sevar (ma:2c4) 

S.var (m11:24) 

SeVvar (ix:24) 

S.evar (ylo324) 

S.evar (yla3z2e4) 

S$evar (ya:24) 

s.fcons (b2,0,168,114,64324) 
s.fcons (61,0,170,114,649:24) 
s.fcons (a2,0,192,164,64224) 
S.1CONns (a1,0,22,228,642:24) 

te. generated forsfile RRAEKKKKKKREKEKEKKKEKKKERKE 
S$eProc Chile? :) 
Seassignconsloneout,0:1,1) 

s.fmul (m2eb,me3,022:224,24,24) 
sefmul (mearmeerbel:24,24,24) 
s.fsub (mearmear,meb:24,24,24) 
s-fsub (melepyle,meas24, 24,24) 
s.fmu] (ybp,ace2e2,m233:24,24,24) 


che 





0000 00 000 0 0 0 0 0 D0 0Z 0 000 


S.tmut 

se fadd 
sefsubd 

S. fx 
Seanaout 
S.fassiaqn 
se-fassign 
seexitoroc 
Sevar 
sevar 
Sevar 
sevar 
Sevar 
Sevar 
sevar 
Sevar 
sefcons 
$.fcons 
s.fcons 
S.fcons 
s.eend 


C(yhracls,m2e2e:24,24,24) 
C(yheyheybs2e4,24,24) 
C(yemeleyhs24,24,24) 
Ciye,y:8) 
Ciy,710,1038) 
(m23,m22324,24) 
(m2ee,me2l13:24,24) 
(fil2s,oneout::) 
(nmeb324) 

mie sted) 

(meazec4) 

(mee:24) 

(yb:24) 

(yh324) 

CD) 

Civ 6) 
(522,0,168,56,64:24) 
(o21,0,36,205,64:24) 
(ace, 0,8,49,64:24) 
(a21,0,4,198,1923:24) 
Gs: 
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APPENDIX B 


DATA, PARALLEL REALIZATION 


This is the input and output data for the parallel 
filter problem developed in chapter five. The timing 
constraints were increased to ninety milliseconds to permit 
successful generation of a realization. The design program 
attempted to generate a dual processor realization for the 
original specification, but a system error prevented 


Sempreetion of the output. 


A. INPUT DATA 
i C€odL Listing 


"This is a parallel realization" 

weaeene fourth order rate filter" 
meleepresented in the article" 

"by Nagle and Nelson in the" 

Beepriary 1931 issued of IEEE COMPUTER." 


meeNlirFiCATION: 
beciamer: M. R. Heilstedt 
Date: 7 April 1983 
Baeqeet: Sample Parallel Filter Problem 


ENVIRONMENT 5 
imme so LL. SDR 1 ,TIL; 
Ouepeessr,o,..L5;CONV,1,TIL; 
Pyagegmeri es y).38:Ml115,8;M12,3;M13.,8, 


Q1A,243Q2A,243;Q1B, 24;Q2B,24;M21,83M22,83M23,8; 
GO ,1;VAR1,1;VAR2,1;CTR1,1;CTR2,1; 


CONTINGENCY LIST: 
"Tnitiate sequence every 60 milliseconds: 
Bore oOms do SETFLG 
"Sample Signal once each period" 
Every, FOms do SAMPLE 


131 6 





"Check for data available after sample converted" 
When DATA:60ms do READ 

"Set flags for parallel operations when ready" 
When LOOP:60ms do FORK 

Seiest Dadrallel filter function" 

When C1:60ms do FILL 

"Second parallel filter function" 

When C2:60ms do FIL2 

"Add results of each function and issue sum" 

Wnen C3:60ms do PLUS 


EROCS DURES; 
Contingency DATA 
Binary DATA,1; 
Sense (DR); 
ie OR=0 then DATA:=1; 
Exit DATA 


Someamaency LOOP 
Banary LOOP,1; 
ihe GO=1 then LOOP:=1; 
GO=0; 
Eee t LOOP 


Pemeimgzency Cl 
pamaray Cl. ls; 
Gav ARi=l then Cl:=1; 
VAR1:=0; 

Sage (eal 


Gomeangency C2 
Beemary C251; 
ioeevek2=i then C2:=sl; 
VAR2:=0; 

eee C2 


uncer .on C3 

Banary C3,1; 

ecu reatr as cand, CER2,g¢.0 then €3:=1; 
Exe C3 


Mas SAMPLE 
CONV:=9; 
Tssue (CONV) ; 
eonVy:=1; 
Tssue (CONV); 

meee on MP LE 





Task READ 
DATA: =03 
Sense(X) ; 
‘CO same 

Exit FORK 


Task FILL 
Cr: =0; 
Pees = X+(15823°M12)-(0.896"°M13); 
Q1B:=€1.264*M11)-(1.748*M12); 
M13;3:=M1235 
Meez 2 =Mi ls 
Oki s=ClTRitl; 

eee ee LL) 


mask FIL2 
C2L=0;5 
Bee =X+ (1. 205°M22)-(0., 444" 23): 
Oye: =(1.673*°M21)=(1.569"M22); 
wees > =M22;3 
Me2:=M21; 
Ome 22 =CTR2+1; 

exit) fiL2 


asic PLUS 
@2--=0'; 
Q1A:=01B; 
¥:=1.0-(8.0*Q1A)-(4.0*02A); 
ieee (Y): 
ers =CiRi-1; 
emp: =CTR2-1: 

exit PLUS 


Blrelts. 
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3. 


Pemecive Listing 


t. generated for:system KKK KKK KKHKKKKKE KKK 
Semain ey 

s.estart (ses 

t. generated forteach KIKI KKK KEK KKK KKK KKK KR 
S.eevery Ceach::) 

Sevar (each:8) 

te generated fors:data RRR KKK KAKKK KEKE KK 
$.OProc (datas: 


s.fixedwait (25:3: 
s.-sensecond (dr:1,1) 


$.eeQ (dte,dr,dc00128,1,1) 
Sejmof ()t2,9d100138) 
Seassigncons(data,1:1,1) 

Se oc (210013: 

Seexitoroc (data,0::) 

Sc Ons (9¢001,031) 

Sevar (Jee: 3) 

S$.var (datas!) 

S$evar Gees tl) 

t. generated for:looo ae ee ee ee ee ee ee ee ee ee ee 
Seproc (loons: 

$.eQq (dt 3,g0,9c002:8,1,1) 
Sejmof (dt 5,010023 8) 


Seassigncons(looo,i:1,!1) 
Seassigncons(go,0:1,1) 


Seloc (d10022:3 

Seexitproc (loop,0:: 
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ihe following errors in the design system were detected 


Bemeorrected during the course of completing this thesis. 


A. DESIGN PROGRAM 

The call to subroutine GETNM for subroutine SYMVAL was 
Corrected. A variable passed in the call, JJ, was 
incorrectly passed with a value of zero. This was changed 
iemeone prior to the call to GETNM. 

Various errors in branch statements existed throughout 
the design program listing and appeared to be the result of 
either an incorrect transfer of the original program to 
magnetic tape or errors in reading the magnetic tape by the 
VAX 11/780. Errors of this nature were detected in sub- 
routines ARGIF, SEDUL, and TMANAL. These have been 
@Semrected, but validation of the remainder of the program 
must be accomplished. 

Subroutine OUTCOM writes the power consumption of the 
realization and its chip count to the terminal screen. This 
was observered to occur twice at the end of each execution. 
Mites subreucine was corrected so that only one such statement 
is now displayed for the user. The destination of the 
duplicate statement was changed to the end of the software 


listing generated for each successful realization. 
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Subroutine SEDUL calls a second subroutine, INSERT, from 
two possible locations. The first of these used an integer 
variable as a parameter for a variable defined within INSERT 
as real. This was observered to prevent successful 
Pe@eration Of realizations for test input. The variable R2Z 
was added as a parameter in the calling routine. It is 
equated to the value of the integer array element 
reek (istop-l,stop) prior to the call to INSERT. 

INSERT sets the values of the array ORDER as part of its 
execution. The type of ORDER is declared to be integer. 
However, the values assigned to it by INSERT are obtained 
from real variables. This caused an integer overflow error 
in the program which resulted in termination of the design 
attempt. The problem was corrected by fixing the values of 
the real variables when equated to the elements of array 


CEE. 


B. REALIZATION LIBRARY 

Entry h.dac was missing the qualifiers on all skip 
statements. These were determined and added. 

Qualifiers for skip statements in EPROM hardware 
Peemitive were set to zero. This caused infinite generation 
of real only memory. The qualifiers were corrected to their 


jpmeper values. 
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Numerous errors attributable to improper transfer of the 
program from magnetic tape to the VAX were detected and 
Memeected. Validation of the complete listing must be 
performed. 

The values of the qualifiers for the SKIP statements in 
primitives h.adce and h.bufframp were incorrect. These were 


corrected to their proper values. 


Se eNCORRECTED ERRORS 

Two errors are known to exist in the main vrogram but 
their source could not be determined. Both of these prevent 
Meeper program execution. 

The first error occurs when a dual processor realization 
memeoenerated. A system error is produced when the program 
attempts to access the file of monitor primitives. 

ihe second error consists of a failure to properly 
evaluate the value of parameters enclosed in the symbols '< 
>', This prevents the proper utilization of library 
primitives which use parameters as qualifiers to logical 


Seaeemencs. 
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